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The Sunyaev-Zel'dovich (SZ) effect from clusters of galaxies is a dominant source 
of secondary cosmic microwave background (CMB) anisotropy in the low-redshift 
universe. We present analytic predictions for the CMB power spectrum from mas- 
sive halos arising from the SZ effect. Since halos are discrete, the power spectrum 
consists of a Poisson and a correlation term. The latter is always smaller than 
the former, which is entirely dominated by nearby bright massive halos, i.e., by 
rich clusters. In practice however, those bright clusters are easy to indentify and 
can thus be subtracted from the map. After this subtraction, the correlation term 
dominates degree-scale fluctuations over the Poisson term, as the main contribu- 
tion to the correlation term comes from distant clusters. We compare the signal 
of the correlation term to the expected sensitivity for the Planck experiment for 
the SZ effect, and find that the correlation term is detectable. Since the degree 
scale spectrum is quite insensitive to the highly uncertain core structures of halos, 
our predictions are robust on these scales. Measuring the correlation term on de- 
gree scales thus cleanly probes the clustering of distant halos. This has not been 
measured yet, mainly because optical and X-ray surveys are not sufficiently sensi- 
tive to include such distant clusters and groups. Our analytic predictions are also 
compared to adiabatic hydrodynamic simulations. The agreement is remarkably 
good, down to ten arcminutes scales, indicating that our predictions are robust for 
the Planck experiment. Below ten arcminute scales, where the details of the core 
structure dominates the power spectrum, our analytic and simulated predictions 
might fail. In the near future, interferometer and bolometer array experiments will 
measure the SZ power spectrum down to arcminutes scales, and yield new insight 
into the physics of the intrahalo medium. 



1 Analytic halo approach to the power spectrum 

We construct an analytic model of the power spectrum of the cosmic microwave 
background (CMB) anisotropy arising from halos through the Sunyaev-Zel'dovich 
(SZ) effect. Since halos, which include clusters of galaxies, groups, and galaxies, 
are discrete objects, their angular power spectrum Ci consists of a Poisson term 
c[ F ^ and a correlation term c[ C ^ that arises from their gravitational clustering^: 

Ci = c\ p) + cf ] . (1) 



szfluctuation'mg9'long: submitted to World Scientific on February 1, 2008 1 




Figure 1. Power spectrum of SZ anisotropics from massive halos (2 X 10 12 h' 1 Mq < M < 
2 X 10 15 h~ 1 Mq) in units of fiK in the Rayleigh- Jeans limit. Thick solid lines show analytic 
predictions for the Poisson and the correlation terms, while dotted lines show the predictions after 
subtracting the ROSAT X-ray flux-limited samples. The number density of subtracted clusters is 
0.8 deg" 2 . Filled and open circles show the power spectrum measured from hydrodynamic simu- 
lations. Filled circles correspond to the dynamic range that is reliably resolved in the simulation. 



The following describes our prescription for computing the SZ power spectrum using 
the halo approach. The key ingredients in our prescription are (1) the halo mass 
function dn(M,z)/dM, (2) the halo bias parameter b(M, z), and (3) the halo, 
pressure profile y(9). (1) and (2) are computed from the Press-Schechter theoryti 
while (3) is calculated assuming that halos are virialized and isothermal. The 
halo approach has also, been used independently to predict the non-linear dark 
matter power spectrumQ and bispectrumtl. In this case, (3) is replaced by the dark 
matter halo profile. These predictions are successful in reproducing the results of 
N-body simulations. Thus, as far as dark matter is concerned, our prescription 
is adequate to describe the halo-halo correlation. However, once gas dynamics 
is included, the degree of uncertainty in (3) increases. We will see below that 
our analytiapiediction nevertheless agrees well with state-of-the-art hydrodynamic 
simulationsBBBEl (see figure g). 
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1.1 Mass function 



One of the key ingredients in our analytic prediction is the halo mass function 
dn(M,z)/dM, which we compute from the Press-Schechter formalismtx This for= 
malism, which plays a central role in our method, was used by Cole and Kaisena 
to calculate the rms SZ fluctuation from halos, including both the Poisson and the 
correlation terms. Makino and SutcO calculated the Poisson term including the 
effect of halo profiles, assuming the f3 profile. Thanks to the recent remarkable 
developments of the microwave experiments, it is now feasible to measure not only 
the rms fluctuation but also the full CMB two point function, or its harmonic trans- 
form, the angular power spectrum Cj—very accurately down to arcminutes angular 
scales. Atrio-Barandela and Mucketo calculated C; from the Poisson term of SZ 
fluctuations using the Press-Schechter mass function and the j3 profile. 

1.2 Bias parameter 

We predict the SZ power spectrum C'i that includes the correlation termi. For this 
purpose, we need the halo-halo correlation function and its evolution with redshift. 
In the framework of the Press-Schechter theory, these can be calculated fronx the 
underlying dark matter correlation multiplied by the bias parameter b(M, z)M. In 
Fourier space, this can be formulated, to the leading order as 

P hh (k, M l7 M 2 , z) ~ b{M u z)b(M 2 , z)P dm (k, z), (2) 

where Phh(k) and Pd m (k) are the halo-halo and the dark matter power spectrum, 
respectively. 

1.3 Halo profile 

In predicting the power spectrum of halos, we need the halo profile which determines 
the small scale power. Since the SZ effect is produced from thermal gas pressure, we 
need to model both the gas density profile and the temperature profile. Firstly, we 
assume that halos are virialized objects with an isothermal temperature distribu- 
tion. Then, we use observationally determined halo gas density profiles, namely the 
/3-profile with /3 = 2/3 and a core radius of 0.15 h' 1 Mpc at M = 10 15 h' 1 M Q smd 
z = 0. We let the core evolve with time according to the self-similar evolutionEj. 
Note that the evolution of the-jcore is the most uncertain parameter in our model. 
Atrio-Barandela and MiicketJi-J employed a different evolution model, and found 
that the arcminute-scales C;'s are sensitive to the core evolution model. We con- 
firmed their result, and also found that the larger angular scale spectrum is fairly 
insensitive to the core model. Thus, as long as angular scales larger than about ten 
arcminutes are considered, uncertainties in core evolution model do not affect the 
prediction. 

2 Power spectrum of the Sunyaev Zel'dovich effect 



szfluctuation'mg9'long: submitted to World Scientific on February 1, 2008 3 



2.1 Poisson and correlation terms 



Using the harmonic transform of the /3 profile for the SZ surface brightness distri- 
bution, yi(M, z), we get the following analytic expressions for the angular power 
spectrum: 

CliP) =^l dZ tl dMdj ^^ M ^^ ® 



dM^^b(M,z) yi (M,z) 



(4) 



where P and C denote the Poisson and the correlation terms, respectively. j v is the 
spectral function of the Sunyaev-Zel'dovich effect and equals —2 in the Rayleigh- 
Jeans regimaij. The range of the mass integral determines what objects are con- 
sidered. We choose 2 x 10 12 h^ 1 M© < M < 2 X 10 15 h^ 1 M , corresponding ap- 
proximately to the resolution and box-size limits of our hydrodynamic simulation!]. 
The redshift integral up to z = 5 is found to be sufficient for convergence. Our 
choice of cosmological parameter is fio — 0.37, Ao = 0.63, fib — 0.049, h = 0.7, and 
cr 8 = 0.8. 

Figure [j] shows the numerical results from equations (|J) and (f|). Firstly, ' 
(p) (p) 

is always smaller than C\ . C\ has a turn over at around I ~ 5000, which 
corresponds to the typical angular size of cope.xadii. Since the spectrum for I > 
2000 is very sensitive to the halo core rnodeliSu, it is highly uncertain. We thus 
mainly focus on larger angular scales (I < 2000). Note however, that the small- 
scale spectrum is of great interest, as it potentially probes the intracluster gas 
pressure state. Although Planck cannot probe the spectrum beyond I ~ 2000, 
forthcoming interferometer and bolometer array experiments will certainly measure 
SZ fluctuations on these scales. 



2.2 Subtracting nearby bright clusters 

Can we measure the correlation term? Yes! Since nearby massive clusters (such as 
the Coma cluster) dominate the Poisson term on large angular scales, the Poisson 
shot noise will be substantially reduced if bright clusters are excised from the map. 
Figure shows the redshift distribution of C; at 3 different angular scales, I = 100, 

1000, and 3000. Nearby halos in z = 0.01 — 0.1 dominate the Poisson C loc ', while 
(c) 

the correlation C 100 come from z = 0.1 — 1. This clear separation in redshift space 
allows us to remove the Poisson shot noise, while preserving the correlation term. 

One way of doing this is to subtract a homogeneous sample of local clusters 
from the SZ map. Figure [l] shows the SZ power spectrum before (solid lines) and 
after (dashed lines) subtracting an X-ray selected cluster sample with the flux limit 
> 10~ 13 erg cm~ 2 s _1 in the ROSAT (0.5 - 2 keV) band. The Poisson spectrum is 
greatly reduced, while the correlation spectrum is less affected by the subtraction. 
The correlation term actually dominates over the Poisson term for I < 200. The 
number density of subtracted clusters is 0.8 deg -2 . For Jiiis figure, the X-ray 
flux was calculated by the method of Kitayama and SutcO. Alternatively, one 
can identify bright SZ clusters in the CMB map itself, and then subtract thenfl. 
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Figure 2. Redshift distribution of the SZ power spectrum in log z-interval, at I = 100, 1000, and 
3000. Amplitudes are normalized to unity. 

For example, 0.9 deg -2 clusters can be removed for ones brighter than 50 mJy at 
350 GHz. 

2.3 Can we detect the correlation term? 

Even if the correlation term overcomes the Poisson term, the primary CMB fluctu- 
ations may prevent us from measuring the correlation term on degree scales. Can 
we really measure the correlation term? Yes! Thanks to the unique shape of the 
frequency spectrum j v of the SZ effectll3, we can subtract the primary CMB, and 
hopefully other foregrounds from the map. In other_svords, we expect to be able to 
derive a cleaned SZ map. Cooray, Hu and TegmarklS demonstrated this technique: 
using the frequency separation, they obtained the expected noise level for the SZ 
power spectrum after subtracting the primary CMB and other modeled sources of 
foreground. Figure || compares the signal of the correlation term to the resulting 
SZ noise power spectrum, binned into bands with Al = 50. The result is so en- 
couraging that we can hope to measure the contribution from the correlation term 
fairly accurately at I = 30 — 200. Since the degree-scale spectrum is insensitive to 
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Figure 3. Comparison of the predicted signal from the correlation term to the expected sensitivity 
of the Planck experiment (filled triangles) to the Sunyaev-Zel'dovich effect. The noise power 
spectrum has been binned into bands with AZ = 50. The thick solid line shows the total spectrum, 
while the dashed lines show the Poisson and the correlation terms separately. The ROSAT X-ray 
flux-limited sample was removed. The number density of subtracted clusters is 0.8 deg - 2 . Note 
that the plotted sensitivity is the residual noise level after subtracting the primary CMB signal 
from the map, on the basis of a multi-frequency analysis. 



the highly uncertain core structures of halos, our predictions are robust on these 
scales. Thus, measuring the correlation term cleanly probes the clustering of dis- 
tant halos. Such a measurement has not been achieved yet, mainly because of the 
lack of sensitivity of optical and X-ray surveys to such distant clusters and groups. 
The fact that the SZ effect is very sensitive to high redshift, is a consequence of the 
known fact that the SZ surface brightness is independent of z. 



3 Comparison to the moving-mesh hydrodynamic simulation 

We carried out adiabatic hydrodxnamic simulations! using the moving-mesh hy- 
drodynamic code written by PenO. The code captures the advantages of both the 
Lagrangian particle based and the Eulerian grid based codes by allowing the grid 
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Figure 4. Snap shot of ^/-parameter in the moving-mesh hydrodynamic simulation at z = 0. The 
box size is L = 100 h^ 1 Mpc, and the number of cells is N = 128 3 . The effective spatial 
resolution at z = is ~ 80 h -1 kpc. Our choice of cosmological parameters is Qo = 0.37, 
Ao = 0.63, f2(, = 0.049, h = 0.7, and as = 0.8. The angular power spectrum of the SZ effect from 
the simulation is plotted in figure BJ. For the simulation, the mean density- weighted temperature 
at present is 0.25 keV, and the mean y-parameter is 1.7 X 10 -6 , in good agreement with the 
Press-Schechter predictions. 



meshes to deform along potential flow lines. This strategy makes it possible to in- 
crease the resolution by twenty fold over previous Cartesian grid Eulerian schemes, 
at a low computational cost. 

Figure ^ shows a snap shot of the y-parameter in the simulation at z = 0. The 
box size is L = 100 h^ 1 Mpc, and the number of cells is N — 128 3 . Since the 
linear compression factor is about 10, the effective spatial resolution at z — is 
~ 80 hr 1 kpc. Our cosmological parameters are fio = 0.37, Ao = 0.63, = 
0.049, h = 0.7, and erg = 0.8. The mass range resolved in this simulation is 
2 x 10 12 h- 1 Ma < M < 2 x 10 15 Z^ 1 M . 
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Figure [l] compares the angular power spectrum of the SZ effect from the sim- 
ulation (filled and open circles) to the analytic predictions (solid lines). The filled 
circles highlight the /-range, 200 < I < 2000, which is reliably described in the 
simulation. The agreement between the simulation and the analytic prediction is 
impressively good. Likewise, the mean density weighted temperature at present 
is 0.25 keV, and the mean y-parameter is 1.7 x 10 -6 , in good agreement with the 
Press-Schechter predictions. Thus, our halo approach to the SZ effect works as long 
as relatively massive halos resolved in this simulation are considered. We conclude 
that the SZ power spectrum in 200 < I < 2000 is dominated by these massive halos. 

4 Discussion 

Our analytic predictions agree well with the adiabatic hydrodynamic simulations; 
however, our halo gas pressure model disagrees somewhat with observations. This 
disagreement affect our analytic predictions and simulated results qualitatively on 
smaller angular scales, say, I > 2000. For example, the self-similar model appears 
to have difficulties in explaining the X-ray luminosity-temperature relatione^, the 
mass-temperature relationEU, and the central entropy-temperature relationE3. The 
departure from the adiabaticity due to so-called preheating is thought to be a 
solution to these problems. If this is so, the adiabatic simulation would not describe 
the intrahalo gas state accurately. The resolution of these disagreements seems to 
require larger core radii than those predicted by the self-similar model. This would 
tend— to suppress the small scale power spectrum, as compared to the adiabatic 
casetil. Including those effects in analytic predictions is rather challenging, but 
is important. Interferometer and bolometer array CMB experiments are expected 
to measure the SZ power spectrum accurately down to arcminutes scales, in the 
near future. These experiments will thus certainly further our understanding of the 
intracluster medium. 
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